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Abstract:- In this paper, Schnackenberg model has been solved numerically for finding an approximate solution 
by Finite difference method and Adomain decomposition method. Example showed that ADM more accurate 
than FDM and more efficient for this kind of problems as shown in tables (1,2) and figures (1-4). 
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I. INTRODUCTION 

Many physical, chemical and engineering problems mathematically can be modeled in the form of 
system of partial differential equations or system of ordinary differential equations. Finding the exact solution 
for the above problems which involve partial differential equations is difficult in some cases. Here we have to 
find the numerical solution of these problems using computers which came into existence [11]. 
For nonlinear partial differential equations, however, the linear superposition principle cannot be applied to 
generate a new solution. So, because most solution methods for linear equations cannot be applied to nonlinear 
equations, there is no general method of finding analytical solutions of nonlinear partial differential equation, 
and numerical techniques are usually required for their solution[6]. 

1.1 Mathematical Model: 

A general class of nonlinear-diffusion system is in the form 

— = t/jAw + fljM + b x v + f(u,v) + gi(x) 

I 

— = d 2 Av + a 2 u + b 2 v - f(u,v) + g 2 {x) 
dt 



With homogenous dirchlet or neumann boundary condition on a bounded domain £2 , n<3, with locally 
lipschitz continuous boundary. It is well known that reaction and diffusion of chemical or biochemical species 
can produce a variety of spatial patterns. This class of reaction diffusion systems includes some significant 
pattern formation equations arising from the modeling of kinetics of chemical or biochemical reactions and from 
the biological pattern formation theory. In this group, the following four systems are typically important and 
serve as mathematical models in physical chemistry and in biology: 

• Brusselator model: 

2 

ay = -(b + 1), by = 0, «2 = ^> ^2 = 0' f = u v > 8\ = a > §2 = 0 
where a and b are positive constants. 

• Gray-Scott model: 

2 

a x = -(F + k), \ = 0, a 2 = 0, b 2 =-F, f = u v, gy = 0, g 2 = F 

where F and k are positive constants. 

• Glycolysis model: 

2 

«1 = -1) by = k, a 2 = 0, b 2 = -k, f = u v, = p, g 2 = 8 

where k, p and 5 are positive constants. 

• Schnackenberg model: 

2 

a\ = ~k, by =a 2 =b 2 = 0, / = u v, gy = a, g 2 = b 
Where k, a and b are positive constants [14]. 

Then one obtains the following system of two nonlinearly coupled reaction-diffusion equations, 
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OH , , 2 

— = a, Au - ku + u v + a 



dt 



(2) 



OV . . 2 , 

— = d~,Av — u v + b 



dt 



With initial and boundary conditions: 

u(t,x) = v(t,x) = 0, t > 0, x £ a/2 ] 
u(0,x) = u 0 (x), v(0, x) — v 0 (x) x £ dH) 
And with Neumann boundary conditions: 

— = — = 0 at x = 0 and x=L (4) 

8x dx 

Where d^,d2 ,k,a,b and L are positive constants [14]. 

Reaction-diffusion (RD) systems arise frequently in the study of chemical and biological phenomena 
and are naturally modeled by parabolic partial differential equations (PDEs). The dynamics of RD systems has 
been the subject of intense research activity over the past decades. The reason is that RD system exhibit very 
rich dynamic behavior including periodic and quasi -periodic solutions [4, 13]. 

Various orders are self -organized far from the chemical equilibrium. The theoretical procedures and 
notions to describe the dynamics of patterns formation have been developed for the last three decades [10]. 
Attempts have also been made to understand morphological orders in biology [5]. Clarification of the 
mechanisms of the formation of orders and the relationship among them has been one of the fundamental 
problems in non-equilibrium statistical physics [9]. 

Various finite difference algorithms or schemes have been presented for the solution of hyperbolic - 
parabolic problem or its simpler derivatives, such as the classical diffusion equation. It is well-known that many 
of these schemes are partially unsatisfactory due to the formation of oscillations and numerical diffusion within 
the solution [12]. 

Solution by the finite difference method, although more general, will involve stability and convergence 
problems, may require special handling of boundary conditions, and may require large computer storage and 
execution time. The problem of numerical dispersion for finite difference solutions is also difficult to overcome 



Adomian decomposition has been applied to solve many functional equations so far. In this article, we 
have used this method to solve the heat equation, which governs on numerous scientific and engineering 
experimentations. Some special cases of the equation are solved as examples to illustrate ability and reliability 
of the method. Restrictions on applying Adomian decomposition method for these equations are discussed [3]. 
The decomposition method can be an inactive procedure for analytical solution of a wide class of dynamical 
systems without linearization or weak nonlinearity assumptions, closure approximations, perturbation theory, or 
restrictive assumptions on stochasticitiy [1]. 



2.1 Finite Difference approximations 

The finite difference Scheme, generally reduces a linear, nonlinear partial differential equations into 
system of linear, nonlinear equations and various methods were developed to find the numerical solution and 

acceleration the convergence [8]. Assume that the rectangle R — {(x,t) : 0 < x < a, 0 < t < b} is subdivided 

into n-1 by m-1 rectangle with sides Ax = h and At =k . Start at the bottom row, where t = t\ = 0, and the 

solution is u(xf, fj) = f(xf) . The grid spacing is uniform in every row: x i+1 — x t + h and (x^ — x t — h), and 

it is uniform in every column: t ;+1 = t ; +feand(t ; _ 1 = t ; — k). And use the approximation U^- for u(x t ,t .) 



[7]. 



II. MATERIALS AND METHODS 



to obtain [8]. 



du 




(5) 
(6) 



dt 

dv 



dt 




u _ u i-l,j- 2u i,j+ u i+l,j 
■ 2 ~ (Ax) 2 
v_ _ v i-l,j- 2v i,j+ v i+l,j 
2 _ (Ax) 2 



(7) 



(8) 



Substitute (5) - (8) in the Schnackenberg model (1) to get 
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u i,j+l ~~ u iJ 
At 

At 



u i,j 



'/•/•I '<,/ 



•2«/,j 



+ u 



(Ax)* 
v i-l,j ~ 2v i,j +v i+l,j 



{AxY 



- hij j + (w ; - j ) v ( - j + fl 



d^At 



(Ax) 
(A*) 2 



2~[";'-l,;' ~ 2u iJ +u i+l,j]~ Atku i,j + At ( u ij) v iJ + 



Ata 



^ v i—l,j ~ 2v iJ + v i+\,j 1 " A/ K, j ) v i, j 



+ Atb 



Let q 



(Ax) z 



and 



(A*) 2 



then 



Kfj+l = «i J + A u i-\,j ~ 2r l"iJ + r l u i+\i ~ Atku i,j + At ( u i,j ) v iJ + Ata 
v iJ+l = v iJ + r 2 v £-l,j " zr 2 v ;,; + r 2 v i+lJ " " 



v ; ;+roV;_i ; -2r 2 Vj-j +r 2 Vj- +1 j -Af(«jj) Vjj+Atb 

2 

u i,j+\ = [ J ~ 2r l ~ A *]";,y + 1 K'-IJ + M i+l,y ] + Af ("; j' ) v ;',y + 



Ata 



v i,j+l = P " 2r 2 l v i,j + r 2 f v /-l , j + v i+l J ] " Af (";j' ) v r J + Atb 

ve 

0 , so u t _i j — and Uij — u 2 j 



(9) 



2 Ax 



= U; 



From the boundary conditions (4) we have 

du _ u i j-u i _ 1 j 
dx ~ 2Ax 

And 

du _ u i+1 j-Ujj 
dx ~ 2Ax 

And also for v 

3v_ _ Vtj-Vj-ij _ 
dx ~ 

And 

Si; _ 
dx 



I ! SO u i + lj 
0 , so = Vi j 

0 , so v l+1J = j 



QYld ^-\\ j — ^10 /' 



and Vij = v 2 j 
and v xl j — v W j 



Unix) 



And from the initial condition: 

"1,1 = "2,1 = "3,1 = "4,1 = "5,1 = "6,1 = "7,1 = "8,1 = "9,1 = "10,1 = "l 1,1 = "0 ( 

v l,l = "2,1 = v 3,l = v 4,l = v 5,l = v 6,l = v 7,l = v 8,l = v 9,l = v 10,l = v ll,l = v oW 

The result equation (8) is the finite difference method for the Schnackenber model. 

2.2 Adomian Decomposition Method 

Nonlinear differential equations are usually arising from mathematical modeling of many frontier 
physical systems. In most cases, analytic solutions of these differential equations are very difficult to achieve. 
Common analytic procedures linearize the system or assume the nonlinearities are relatively insignificant. Such 
procedures change the actual problem to make it tractable by the conventional methods. This changes, some 
times seriously, the solution. The above drawbacks of linearization and numerical methods arise the need to 
search for alternative techniques to solve the nonlinear differential equations, namely, the analytic solution 
methods, such as the perturbation method, the iteration variational method [11-14] and the Adomian 
decomposition method. 

The decomposition method was first introduced by Adomian since the beginning of the 1980s. The 
Adomian decomposition method (ADM) is used to solve a wide range of physical problems. This method 
provides a direct scheme for solving linear and nonlinear deterministic and stochastic equations without the need 
for linearization this yields convergent series solutions rapidly. An advantage of this method is that, it can 
provide analytical approximation or an approximated solution to a rather wide class of nonlinear (and stochastic) 
equations without linearization, perturbation, closure approximation, or discretization methods. Unlike the 
common methods, i.e., weak nonlinearity and small perturbation which change the physics of the problem due 
to simplification, ADM gives the approximated solution of the problem without any simplification. Thus, its 
results are more realistic [1,2]. 



We define the operator L t 



-1 



CI 



: \{.)dt 

0 



and L Y 



then 



www.ijeijournal.com 



Page 173 



Some numerical methods for Schnackenberg model 



system (1) can be written as: 

2 

L t u = d^L-^u - ku + u v + a 



L t v = d 2 L xx v — u v + b 



By applying the inverse of operator L t on (10,1 1) we get: 

u(t, x) = u(0, x) + d x L~ l (L xx u) - kL, l u + L' 1 (u 2 v) + L^a 

v(t,x) = v(0, x) + d 2 L l l (L xx v) - L, 1 (u 2 v) + L~ l b 

By initial conditions in system (3) then (10,1 1) can be written as: 

u(t,x) = uq(x) + d^Lf ^(L^u) -kL t + L t (u v) + L t 

-1 -1 2 -1 

v(t,x) = Vq(x) + d 2 L( (A*x v ) _ ^f ( M v ) + ^f b 

By using Adomian decomposition method : 



(10) 
(11) 



u(t,x)= X Uf,(t,x) and v(f,jc)= E v„(/,x) 
n=0 n=0 



00 _j CO _j 00 _| OD 

E u n (t,x) = uq(x) + d^L t {Lfx E u n )-kL t u n + L t X A. n +L t L a 
n=0 n=0 n=0 n=0 



r-1, 



r -l 



r-1, 



I v n (t,x) = v 0 (x) + d 2 L t (L^ Z v„)-L, ^ B n + L t b 
n=0 n=0 n=0 



(12) 



(13) 



Where A n and A M are Adomian polynomials. But A n =B n because both non-linear terms are U V 



Where A n 



1 d 



n\dX n 



i'=0 



1=0 



But here A„ 



1 <r 



F( E A'm,-, E A'vj) 



n\ dX n L (=0 (=0 _U=0 
then by equation (12): 
"0 ="0^ 

u k+l = d 1 Lj 1 (L xx u k )-kLj 1 u k +L 7 lA yt+ L 7 la 
By equation (13): 
v Q =v 0 (x) 

v k+l = d 2 L 7 ( L xx v k )- L i B k +L ~t b 



because non-linear term have two functions u(t,x) and v(t,x) 



k>0 



k>0 



A Q =- 



1 d' 



X=0 



0!rf2°L i=0 l 'i=0 
2 

= F(Uq,Vq) = UqVq =B 

k=0 

"1 = d l L J l ( L xx u 0^~ kL t lu 0 + L l LA Q+ L V a = d lk' ( l xx u q) ~ kh^O + L 7 1 ("O v 0^ fL 7 la 

2 2 
= d^L^UQt - hiQt + UQVQt + at = [d^L^UQ - Iqiq + Uq Vq + a]t 



(14) 



,-1 



-I, 



,-1 



-I, 2 



Le? : t/1 = d^L^UQ - hiQ + uqVq + a 
Mr =Ult 

vi = d 2 Lj 1 (L xx v 0 )- Lj l B Q + Lj { b = d 2 Lj l (Lxx-vq) - Lj 1 (uqvq) + Lj { b 

2 2 
= d 2 L xx v 0< ~ U 0 V 0< + bt = t d 2 L xx v 0 ~ "0 V 0 + W 
2 

Let:V\ = d 2 L xx v Q -uqvq +b 
v. = Vlt 



1 d 



V-dX l l 'i=0 



1 i 1 i 

F( ZZ_X Mj, E 1 V;) 



!=0 



A=0 



[F(kq +^1^0+^)]^^ 



(15) 



(16) 
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d 2 d 2 2 2 

= — K"0 +/l "l) (yQ+Xv x )\ x= Q =— [(uq + 2Xu q u x +1 u x )(v 0 +Avi)]^ =0 
dX dX 

d 2 22 22 32 

= — [uqVq+22mqu^vq+ X «j vq + /IwqVj + 2/1 uqu^v^+X u \ v \\x=0 
dX 

2 2 

= 2«qKj Vq + Wq Vj = 2«q (£/ 1?)Vq + «Q Vlf 

Aj = [2m 0 v 0 [/1 + UQVljt = (17) 
K=l 

«2 - ^i^t ^ (Lxx u l ) ~~ ^ ^"1 + ^1 + = ^l^f ^[Ljx (Ult)] - kL t ^ (Ult) + L t ^ ([2wqvq[/1 + uq Vl]t) 

2 2 2 f 2 

= d l L xx (UV) — -kUl — + [2u 0 v Q Ul + VI] — = [d^ {U\)-MJ\+ IuqVqUI + Kg VI] — + at 

Let:U2 = d^L^ (Ul)-kUl + 2u Q v Q Ul + uqVI 
t 2 

w 9 =t/2 — + at (18) 
z 2 

v 2 = d 2 Lj l (Z^v,) + FLj 1 (1- v,) - L7 1 Sj 

= d 2 L 7 1 [^xc (VI*)] + ^ (1 - VI/) - Lj\[2u Q v 0 Ul + kqVi]?) 

2 2 2 

= (VI) y + F(t - VI — ) - [2u 0 v 0 Ul + uqVI] — 

2 t 2 

v 2 =Ft + [rf 2 L xc (VI) -FV\- 2u q vqU\ - uq VI] — + bt 

Let : V2 = d 2 L xx (VI) - FV1 - 2u Q v Q Ul - u 2 Vl 
t 2 

v ? =Ft + V2 — +bt (19) 
L 2 



1 d 2 
2! dX Z 



2 ■ 2 ■ 

F( £ X Uj , £ X V; ) 
!=0 !=0 



X=0 21 dX 



1 " 2 2 
2"[F(mq + X u 2 ,vq + Xv^ + X v 2 )]^ = q 



1 j2 

1 2 2 2 

-^-[{uq + Xu x + X u 2 ) v 2 )]l=0 



2!^ 
1 rf 2 

^-[[uq +2uq(Xu^+ X u 2 ) + (Xu^+ X u 2 ) ](vq+ Xv^+ X v 2 )]^_q 

2 - dX 

1 d 2 

— — ^-[(uq +2Xuqu^+2X uqu 2 +X uy +2X u^u 2 +X u 2 )(vq+Xv^+X v 2 )]^_q 
i j2 

1 w 2 2 223 42 

y [mqVq + 2/ImqKj Vq + 2/1 UqU 2 Vq + X Mj Vq + 2/1 U^U^Q + X "2 V 0 



i j2 

1 d rr 2 , „ . ,2 , . , ,2 x 2 n , .0 . ,2 



22 3 324 52 

+ Xuqv^+2X uqu^ vj + 2A uqu 2 v^ + X v^ + 2/1 u^u 2 v^+Xu 2 v^ 

^ ^ 4 2 5 6 2 

+ A mqV2+2/1 «qMjV2+2/1 MQM2V2+/I MjV2+2/l wj«2V2+/i. "2 y 2 1/1=0 

1 2 2 2 2 

= — (4«Qi<2VQ + 2«j Vq + 4wqKjVj + 2«qV 2 ) = 2uqU 2 Vq + 2wqMjVj + «QV2 + «J Vq (20) 

2 2 
A 2 =2u Q (U2— + at)v Q + 2mq (Ult + at)(V\t + bt) + u 2 (Ft + V2—) + (Ult + at) 2 v Q 

By eqs (16), (17), (19) and (20) 

2 

A 2 = i/ 0 v 0 ;72f 2 + 2u Q UlVlt 2 + u 2 V2— + (Ul) 2 v Q t 2 + (uQbt2u Q a)t 

A 2 =[M 0 v 0 (/2 + 2M 0 (/lVl + i;<QV'2 + (f/l) 2 v 0 ]f 2 + (u 2 bt2u Q a)t = B 2 (21) 
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k = 2 

«3 = d\L t ^ (L XX U2 ) - kLf ^w 2 + L^a 

By equations (18) and (21) 

2 2 j 

k 3 = rf 1 L7 1 [L^ (1/2 — )] - H^ 1 (t/2— ) + Lj 1 ([u Q v 0 U2 + 2u Q UlVl + - UqV2 

+ (Ul) 2 v Q ]t 2 ) 

3 3 3 

= -d x L xx (U2) t —--kU2 t — + [u 0 v 0 U2 + 2u Q UlVl + - uqV2 + (Ul) 2 v Q ] *— 

1 1 1 2 2 t 3 

u 3 = [- dyh^ (t/2) --W2 + u 0 v q U2 + 2u Q UlVl + ~uqV2 + (Ul) vq ] y 

2 ^ 

+ (a + u^b + 2uqO) h ar) (22) 

v 3 =rf 2 L7 1 (L xc v 2 )-L7 1 S 2 +L7 1 * 
Byeqs(19)and(21) 

f 2 1 
v 3 = rf 2 L7 ! [L^ (V2 — )] - L7 1 ([M 0 v 0 f/2 + 2k 0 (/1V1 + - UqV2 + (Ul) 2 v Q ]t 2 ) + u^b + 2u Q a)t) + bt 

3 3 9 

1 f 1 ? 2 t 9 f 

= - d 2 L xx ( v 2 ) — " IuqVqU2 + lu ()UWl + - uqV2 + (Ul) vq ] — + (u^b + 2u Q a) — + bt 

3 2 

v 3 = [- (V2) - uq V qU2 - 2u Q UlVl - - uqV2 - (Ul) 2 v Q ] '— + (u^b + 2u Q a) + bt (23) 



III. APPLICATION (NUMERICAL EXAMPLE) 

We solved the following example numerically to illustrate efficiency of the presented methods. 

— = diAu-ku + u 2 v + a ,t>0, x £ H 

at 1 

— = d, Av - u z v + b ,t>0, x 6 12 
We the initial conditions 

u(x, 0) = Us + 0.01 sm(xx/L) for 0 <x <L 

v (x, 0) = Vs - 0. 12 sin(nx/L) for 0 <x <L 

w(0, t) = Us, u(L, t) = Us and v (0, t) = Vs, v (L, t) = Vs 

We will take 

d,=d 2 =0.01 , a=b=0.09 ,k=-0.004, Us=0, Vs=l 



IV. FIGURES AND TABLES 
Table 1 Comparison between the FDM and ADM for the values of concentration V. 



t= 1 



t=2 



t=3 



ADM FDM ADM FDM ADM FDM 



0 


1.0001 


1.0000 


1.0001 


1.0000 


1.0035 


1.0000 


0.1 


1.0014 


1.0008 


1.0022 


0.9993 


1.0048 


0.9971 


0.2 


1.0020 


1.0010 


1.0022 


0.9976 


1.0043 


0.9864 


0.3 


1.0021 


1.0019 


1.0012 


0.9946 


1.0021 


0.9787 


0.4 


1.0018 


1.0012 


0.9994 


0.9906 


0.9984 


0.9704 


0.5 


1.0011 


1.0005 


0.9971 


0.9861 


0.9936 


0.9624 


0.6 


1.0003 


0.9996 


0.9945 


0.9817 


0.9881 


0.9553 


0.7 


0.9995 


0.9982 


0.9921 


0.9777 


0.9826 


0.9499 


0.8 


0.9988 


0.9968 


0.9900 


0.9746 


0.9781 


0.9464 


0.9 


0.9984 


0.9954 


0.9887 


0.9726 


0.9750 


0.9452 


1.0 


0.9982 


0.9942 


0.9882 


0.9719 


0.9739 


0.9464 


1.1 


0.9984 


0.9954 


0.9887 


0.9726 


0.9750 


0.9499 
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1 ^1 

1.2 


a nnoo 

0.9988 


A AA/C O 

0.9968 


A AAAA 

0.9900 


a rn a c 

0.9746 


A A"7 O 1 

0.9781 


A AC 

0.9553 


1.3 


0.9995 


a a a o ^ 

0.9982 


0.9921 


0.9777 


A A O T7 

0.9827 


0.9624 


1.4 


1.0003 


0.9996 


0.9945 


0.9817 


0.9881 


0.9704 


1.5 


1.0011 


1.0005 


0.9971 


0.9861 


0.9936 


0.9787 


1.6 


1.0018 


1.0012 


0.9994 


0.9906 


0.9984 


0.9864 


1.7 


1.0021 


1.0019 


1.0012 


0.9946 


1.0021 


0.9927 


1.8 


1.0020 


1.0010 


1.0022 


0.9976 


1.0043 


0.9971 


1.9 


1.0014 


1.0008 


1.0022 


0.9993 


1.0048 


0.9993 


2.0 


1.0001 


1.0000 


1.0010 


1.0000 


1.0035 


1.0000 



Table 2 Comparison between the FDM and ADM for the values of concentration U. 





t=l 




t=2 




t=3 




X 


ADM 


FDM 


ADM 


FDM 


ADM 


FDM 


0 


0.0015 


0 


0.0048 


0 


0.0081 


0 


0.1 


-0.0150 


-0.0232 


-0.0102 


-0.0276 


-0.0108 


-0.0319 


0.2 


-0.0306 


-0.0456 


-0.0241 


-0.0539 


-0.0221 


-0.0623 


0.3 


-0.0450 


-0.0664 


-0.0367 


-0.0783 


-0.0316 


-0.0902 


0.4 


-0.0579 


-0.0851 


-0.0476 


-0.0999 


-0.0390 


-0.1149 


0.5 


-0.0692 


-0.1014 


-0.0568 


0.1185 


-0.0444 


-0.1359 


0.6 


-0.0786 


-0.1151 


-0.0643 


-0.1339 


-0.0480 


-0.1532 


0.7 


-0.0861 


-0.1258 


-0.0700 


-0.1459 


-0.0502 


-0.1666 


0.8 


-0.0915 


-0.1336 


-0.0741 


-0.1545 


-0.0514 


-0.1761 


0.9 


-0.0948 


-0.1383 


-0.0766 


-0.1596 


-0.0520 


-0.1818 


1.0 


-0.0959 


-0.1398 


-0.0776 


-0.1613 


-0.0522 


-0.1837 


1.1 


-0.0949 


-0.1383 


-0.0770 


-0.1596 


-0.0523 


-0.1818 


1.2 


-0.0916 


-0.1336 


-0.0749 


-0.1545 


-0.0519 


-0.1761 


1.3 


-0.0863 


-0.1258 


-0.0710 


-0.1459 


-0.0509 


-0.1666 


1.4 


-0.0788 


-0.1151 


-0.0654 


-0.1339 


-0.0488 


-0.1532 


1.5 


-0.0693 


-0.1014 


-0.0578 


0.1185 


-0.0451 


-0.1359 


1.6 


-0.0580 


-0.0851 


-0.0484 


0.1185 


-0.0395 


-0.1149 


1.7 


-0.0450 


-0.0664 


-0.0372 


-0.0999 


-0.0319 


-0.0902 


1.8 


-0.0306 


-0.0456 


-0.0244 


-0.0783 


-0.0223 


-0.0623 


1.9 


-0.0150 


-0.0232 


-0.0103 


-0.0539 


-0.0108 


-0.0319 


2.0 


0.0015 


-0.0000 


0.0048 


-0.0000 


0.0023 


-0.0000 
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V. CONCLUSION 

The Schnackenberg model solved Numerically using finite difference method and Adomain 
decomposition method, and we found that's finite difference method is earlier that ADM but ADM is more 
accurate than FDM and more efficient as show in tables ( 1-2 ) and figures (1-4 ). 
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